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ABSTRACT 

A triplet of subordinate lines of Mg II exists in the region around the h&k lines. In solar spectra these lines 
are seen mostly in absorption, but in some cases can become emission lines. The aim of this work is to study 
the formation of this triplet, and investigate any diagnostic value they can bring. Using 3D radiative magneto¬ 
hydrodynamic simulations of quiet Sun and flaring flux emergence, we synthesize spectra and investigate how 
spectral features respond to the underlying atmosphere. We find that emission in the lines is rare and is typically 
caused by a steep temperature increase in the lower chromosphere (above 1500 K, with electron densities above 
10 18 m -3 ). In both simulations the lines are sensitive to temperature increases taking place at column masses 
> 5 • 10 -4 gem -2 . Additional information can also be inferred from the peak-to-wing ratio and shape of the 
line profiles. Using observations from NASA’s Interface Region Imaging Spectrograph we find both absorption 
and emission line profiles with similar shapes to the synthetic spectra, which suggests that these lines represent 
a useful diagnostic that complements the Mg II h&k lines. 

Subject headings: Sun: atmosphere — Sun: chromosphere — radiative transfer 


1. INTRODUCTION 

The Mg II h&k resonance lines are among the strongest in 
the solar spectrum. They are formed higher than the widely- 
studied Ca II H & K lines owing to the higher magnesium 
abundance. However, because they sit in the UV spectrum 
they are not observable from the ground and have not been 
routinely observed in the past. To observe them astronomers 
have used a multitude of space platforms, from balloons to 
space missions (e.g. Durand et al. 1949; Bates et al. 1969; 
Staath & Lemaire 1995; Doschek & Feldman 1977; Morrill 
& Korendyke 2008; West et al. 2011). The advent of the 
IRIS mission (De Pontieu et al. 2014) has provided unprece¬ 
dented continuous time series of Mg II spectra (and slit-jaw 
filtergrams) with high spatial, spectral, and temporal resolu¬ 
tion, which no observatory before could do concurrently. This 
wealth of Mg II spectra is making a difference in how the chro¬ 
mosphere is understood and has great potential for the future. 

To understand the complex formation of the h&k lines in 
the Sun, several studies have been undertaken, starting with 
the early work of Dumont (1967), Milkey & Mihalas (1974), 
Ayres & Linsky (1976), Gouttebroze (1977), and Uitenbroek 
(1997). More recent work has focused on understanding their 
formation over a range of solar positions (Avrett et al. 2013), 
their polarization potential (Belluzzi & Trujillo Bueno 2012), 
their formation in prominences (Heinzel et al. 2014), and their 
diagnostic value using 3D radiative magneto-hydrodynamic 
(rMHD) models (Leenaarts et al. 2013a, b; Pereira et al. 2013). 
All of these studies focused on the h&k lines; notably ab¬ 
sent was any detailed investigation on the companion triplet 
of lines between the 3 p 2 P and 3 d 2 D states (an exception is 
Feldman & Doschek 1977, who studied these lines above the 
solar limb), with vacuum wavelengths of 279.160, 279.875, 
and 279.882 nm (henceforth referred to as triplet lines). 

With their lower levels being the upper levels of the h&k 
lines (see Figure 1 of Leenaarts et al. 2013a), the triplet lines 

tiago .pereira @ astro .uio. no 


sit on the wings of the h&k lines (one line on the blue wing 
of the k line, and two overlapping lines located between the 
k and the h line). This transition structure bears some re¬ 
semblance to that of the Ca II atom, whose infrared triplet 
lines (849.8, 854.2, and 866.2 nm) share the upper level of the 
H&K lines. With the Ca II infrared triplet firmly established as 
a chromospheric diagnostic in recent literature (in particular 
the 854.2 nm line, see e.g. Cauzzi et al. 2008, Reardon et al. 
2009, Leenaarts et al. 2009, de la Cruz Rodriguez et al. 2013 
and references therein), a study of the diagnostic potential of 
the Mg II triplet is both timely and relevant. 

The Mg II triplet lines are generally much weaker than their 
famous siblings h&k and appear mostly as absorption lines. 
Nevertheless, in energetic events they become emission lines. 
Given the gap in the literature on these lines, in this work we 
set forth to answer the following questions: in which condi¬ 
tions do these triplet lines form, and what can we learn from 
them? Under what circumstances do they become emission 
lines? The lines have already been extensively observed by 
IRIS (the lines around 279.88 nm are included in virtually all 
IRIS observations), so any insight can be useful for a wide 
range of observations. To understand the formation of the 
lines we employ 3D rMHD models, in a similar way to what 
Leenaarts et al. (2013b) and Pereira et al. (2013) did for the 
h&k lines. 

The outline of this paper is as follows. In §2, we describe 
the simulations used and how the synthetic spectra were cal¬ 
culated, and in §3 we study how the Mg II triplet lines are typi¬ 
cally formed in quiet Sun. In §4 we investigate the conditions 
that lead to emission in the triplet lines, both in a quiet Sun 
simulation and a flaring simulation. In §5 we show some ex¬ 
amples of emission profiles observed with IRIS, and we con¬ 
clude with a discussion in §6. 

2. SYNTHETIC SPECTRA 

To study the formation of the Mg II triplet, we follow the 
approach of Leenaarts et al. (2013b) and Pereira et al. (2013) 
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and make use of 3D rMHD simulations performed with the 
Bifrost code (Gudiksen et al. 201 1). 

Bifrost solves the resistive MHD equations on a staggered 
Cartesian grid. A 24 x 24 x 16.8 Mm 3 region of the solar at¬ 
mosphere was simulated, with a constant horizontal cell spac¬ 
ing of 48 km and non-uniform vertical spacing, extending 
from 2.4 Mm below the average T 500 = 1 height and up to 
14.4 Mm, covering the upper convection zone, photosphere, 
chromosphere, and lower corona. The radial curvature of the 
Sun is neglected. The simulations employed include a de¬ 
tailed radiative transfer treatment including coherent scatter¬ 
ing (Hayek et al. 2010), a recipe for non-local thermodynam¬ 
ical equilibrium (non-LTE) radiative losses from the upper 
chromosphere to the corona (Carlsson & Leenaarts 2012), and 
thermal conduction along magnetic field lines (Gudiksen et al. 
2011). We use two simulation snapshots from different runs: 
an “enhanced-network” quiet Sun simulation and an emerging 
flux simulation with some small flares. The quiet Sun simula¬ 
tion snapshot is the same as used in other papers of this series 
(see Leenaarts et al. 2013b, and references therein for more 
details), which includes non-equilibrium hydrogen ionization 
in the equation of state (Leenaarts et al. 2007). The photo- 
spheric mean unsigned magnetic field strength of the simu¬ 
lation is about 5 mT (50 G), concentrated in two clusters of 
opposite polarity, placed diagonally 8 Mm apart in the hori¬ 
zontal plane. 

We also use a snapshot from the “flaring” simulation of 
Archontis & Hansteen (2014). In this emerging flux simu¬ 
lation a uniform magnetic flux sheet with B y = 336 mT and 
of dimension 24 x 13Mm 2 is injected into a numerical do¬ 
main of 24 x 24 x 17 Mm that contains a weak initial field 
of B < 10 -5 mT. The model has fully developed convection 
and a certain percentage of the injected field emerges into the 
chromosphere and corona in a non uniform manner, leading 
to patchy reconnection, as loops expanding through the photo¬ 
sphere into the upper atmosphere come into contact with each 
other. The reconnection leads to structures resembling small 
flares. We use the snapshot at t = 10000 s in which several 
such flares are present. Lurther details of this simulation setup 
can be found in Archontis & Hansteen (2014). 

The synthetic spectra were calculated using the RH1.5D 
code (Pereira & Uitenbroek 2015), a modification of the RH 
code (Uitenbroek 2001) that solves the non-LTE problem for 
each column in a 3D atmosphere as an independent ID col¬ 
umn. As shown by Leenaarts et al. (2013a) this is a good 
approximation for the Mg II h&k lines, outside the I 13 and 
k 3 cores. To reduce the computational costs, we performed 
the calculations for every other spatial point in the horizon¬ 
tal directions. We find that the effects of partial redistribu¬ 
tion (PRD) in the triplet lines are negligible, therefore we as¬ 
sumed complete redistribution (CRD) for all calculations of 
these lines, while calculating the h&k lines in PRD (using the 
fast angle-dependent approximation of Leenaarts et al. 2012). 

3. FORMATION OF THE MG II TRIPLET IN QUIET SUN 

In the quiet Sun, the Mg II triplet lines are observed mostly 
as absorption lines (see e.g. Morrill & Korendyke 2008). This 
is also the case in our synthetic spectra from the quiet Sun 
simulation: the lines are in emission in « 0.5% of the points. 
The lines are formed at around 0.6 — 1.2 Mm above the height 
where T 500 = L The 279.160 nm line is formed at slightly 
lower heights than the other pair, which are blended and are 
indistinguishable at the spectral resolution of IRIS - they usu¬ 


ally appear as a single, wide absorption feature. Compared 
to the 279.160 nm line, the 279.882 nm line has an oscillator 
strength twice as high; in terms of formation height this trans¬ 
lates to a difference of less than 100 km in most cases. The 
279.875 nm line is the weakest, which makes the blend in this 
region asymmetric, with a centroid shifted toward the red. 

The lines exist in a heavily blended region. We synthesized 
the strongest nearby lines (assuming LTE) and find that the 
contribution of lines other than Mg II is negligible through¬ 
out the line profiles, except in the far wings of the Mg II 
279.160 nm line. 

In quiet Sun conditions, the triplet lines are unremarkable 
compared to other strong lines in the region (see Pereira et al. 
2013). They probe similar atmospheric layers to other strong 
lines in between the h&k lines. The line pair at 279.88 nm is 
problematic to measure because the lines overlap. Its asym¬ 
metric shape complicates the task of measuring velocities or 
widths with standard techniques (e.g. position of centroid or 
line fitting). Strong shifts can also make the lines overlap with 
other nearby lines. In Figure 1 we analyze the formation of 
the 279.88 nm lines in two columns from the quiet Sun sim¬ 
ulation, in the format developed by Carlsson & Stein (1997). 
Case (a) shows a typical quiet Sun profile. The two over¬ 
lapping triplet lines at 279.88 nm are the two humps around 
Av = 0, while the feature at Av ~ 23 kms -1 is a Fe I line. 
Here it can be seen that the source function loosely follows 
the Planck function up to a height of about 0.3 Mm, and then 
continues to drop, causing an absorption line. (For this lo¬ 
cation, the h&k lines decouple from the source function at 
z ~ 0.8 Mm, and still follow the lower part of the tempera¬ 
ture increase until z ~ 1.5 Mm.) Case (b) shows one of the 
rare cases of emission in the quiet Sun simulation. Here the 
source function dips with the temperature minimum around 
z ~ 0.3 Mm, but then follows a temperature increase and 
peaks at around z ^ 0.7 Mm, causing an emission line. 

In Figure 2 we compare the formation region of the two 
strongest triplet lines with that of the Mg II k 3 and k 2 features 
(calculated in the same manner as in Leenaarts et al. 2013b) 
for the quiet Sun simulation. The distributions are computed 
using Kernel Density Estimation (KDE, see Rosenblatt 1956, 
Parzen 1962) using a Gaussian kernel. In the case of k 2 we 
took the average z(T = 1) and column mass of the k 2 v and 
k 2 r features. One can see how the Mg Ilk features and the 
triplet lines cover distinct regions in the range of 0.5 — 3 Mm, 
or 10 -1 — 3 • 10 -7 gem -2 and therefore can complement each 
other as diagnostics of the chromosphere. 

4. EMISSION IN THE MG II TRIPLET 
4.1. Conditions for emission 

Under particular conditions the triplet lines become emis¬ 
sion lines. This seems to happen when there is a rapid increase 
of temperature with height in the region around the tempera¬ 
ture minimum. With such temperature rises, the source func¬ 
tion is still close to the Planck function and follows its rise 
before dropping down in higher layers, causing an emission 
line. In the synthetic spectra, emission in the line cores is 
very rare in the quiet Sun simulation (« 0.5% of the columns) 
but much more common in the flaring simulation (« 90% of 
the columns). The conditions for emission become clear in 
Figure 1, where in panel (b) a steep temperature increase at 
0.3 < z < 0.8 Mm leads the source function to peak and then 
drop, causing an emission line. 

Here we define “emission” as the intensity in the line core 
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Figure 1 . Intensity formation diagram for the two Mg II lines around 279.88 nm, in two columns of the quiet Sun simulation. Case (a) shows a typical column 
(left side), with the lines in absorption, while case (b) shows a column where the lines are in emission (right side). For (a) and (b), the individual panels show the 
quantity specified in the top left corner, where Xv is the absorption coefficient, T v the optical depth, and S v the total source function. The quantities are shown as 
a function of Doppler shift from 279.88 nm and simulation height z. The top left panels (Xv Av ) represent the asymmetry contribution. The bottom right panel 
shows C/, the contribution function to intensity, and is obtained by multiplying the other three panels (to improve its visibility, Q was divided by the maximum 
at each wavelength). A t = 1 curve (cyan) and the vertical velocity (red, positive is upflow) are plotted in each panel, with a v z = 0 line in the first panel for 
reference. In the upper right panel we show also the Planck function (white dotted) and the source function at Av = 0 (white dashed) in temperature units (scale 
at the top). The lower right panel contains also the intensity profile in brightness temperature units (scale on the right), at the resolution of the simulation (solid 
white line) and convolved with the instrumental profile of IRIS (dashed white line, spatially and spectrally convolved as in Pereira et al. 2013). 



Height (Mm) 



Column Mass (g cm 2 ) 

Figure 2. Comparison of quiet Sun formation regions for the Mg II triplet 
lines and the k line. Showing distributions of the height where optical depth 
reaches unity and corresponding column mass, for the line center wave¬ 
lengths of two triplet lines, and for the k 3 and k 2 features of the k line. Curves 
depict the kernel density estimation with a Gaussian kernel (see text), normal¬ 
ized by the maximum value. 

being higher than in the line wings. In the quiet Sun sim¬ 
ulation one finds about 1% of the columns with emission 
“bumps” in the far wings - these correspond to locations 


where the heating takes place at deeper layers and have es¬ 
sentially the same formation mechanism as the profiles with 
emission in the line core. Nevertheless, for the remainder of 
this section we restrict ourselves to the more extreme events 
of line core emission. 

In Figure 3 we show distributions for the ranges of heights, 
column mass, and electron densities where the lines are 
formed. Unlike in Figure 2, where we show the distributions 
for z(t = 1 ), here we show the ranges of typical conditions 
where the lines are formed. The distributions are given for the 
maximum and minimum z (T = 1), the heights where the op¬ 
tical depth reaches unity, and for the values of column mass, 
electron density and temperature at those maximum and min¬ 
imum heights. The bulk of the line is formed in the region be¬ 
tween the maximum and minimum distributions, with the line 
wings formed close to the minimum heights, and the line cen¬ 
ter formed closer to the maximum heights. Reflecting the very 
different density profiles, the lines are formed much higher 
in the emerging flux flaring simulation than in the quiet Sun 
simulation. However, when shown on a column mass scale 
the formation region is very similar. The electron density and 
temperature distributions are more convoluted, and again re¬ 
flect the differences between the simulations. In the flaring 
simulation there is a subset of columns with very cool mid¬ 
chromospheres, and this causes a double peaked distribution 
for the temperature of maximum height of formation. There 
is a clear tendency for the regions in emission to be formed 
at lower column mass densities, higher up in the atmosphere; 
this is true for both the quiet Sun simulation and the flaring 
simulation. We find that the triplet lines are sensitive to col- 
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Figure 3. Distributions of physical conditions for the maximum and mini¬ 
mum height of unity optical depth, for a spectral range of ±30 kms -1 around 
the two strongest Mg II triplet lines. Curves depict the kernel density estima¬ 
tion with a Gaussian kernel, normalized by the maximum value. Results are 
shown for the quiet Sun (QS) simulation for all columns (blue), and for the 
columns with emission in the triplet lines (yellow), as well as for the flaring 
simulation (all columns, black). 

umn masses of 5 • 10~ 4 g cm -2 and higher, in regions with 
electron densities between 10 16 — 10 19 m -3 . 

The emission in the lines is caused by a source function 
increase in the formation region of the line centers. Typi¬ 
cally, this increase in the source function is caused by a large 
temperature increase in the lower chromosphere (> 1500 K). 
Nevertheless, there are many cases when large increases do 
occur, but the source function does not follow the tempera¬ 
ture and the lines are in absorption. While there are many fac¬ 
tors that affect the coupling between the source function and 
the local temperatures, the source function tends to be more 
tightly coupled to the temperature in regions with higher elec¬ 
tron density. In Figure 4 we show a scatter plot and distribu¬ 
tions of AT versus the maximum electron density, for the quiet 



Maximum electron density (m 3 ) 


Figure 4. Physical conditions for emission in the strongest Mg II triplet lines, 
from the quiet Sun simulation. Temperature difference AT in the line forma¬ 
tion region (see text) is plotted against the maximum electron density in the 
region where line core is formed. Showing all simulation locations (small 
blue dots) and locations showing emission in the lines (red circles). The top 
and right diagrams show the distributions (Gaussian kernel density estima¬ 
tion, normalized) for the electron density and AT, respectively. 
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Figure 5. Relation between the ratio of line core to wing intensity (for the 
strongest of the Mg II triplet lines) vs. AT, a temperature difference cor¬ 
responding to Smax — T m in (blue triangles) or T(S max ) — T m - m (red circles), 
where T m j n is the temperature minimum between column masses of 10 to 
0.1 g cm -2 , Smax is the maximum of the source function (in brightness tem¬ 
perature units) between column masses of 0.5 to 10~ 3 g cm -2 , and T(S max ) 
is the temperature at the height where the source function maximum occurs 
(all points for locations in emission). 

Sun simulation. For each simulation column, AT is the differ¬ 
ence between the maximum temperature between the column 
masses of 5 • 10 -4 and 10 -2 gem -2 and the minimum temper¬ 
ature between the column masses of 10 -1 and 10 gem -2 (in 
other words, the temperature difference between the line core 
forming regions and the line wing forming regions). The max¬ 
imum electron density n e is taken between the column masses 
of 5 • 10 -4 to 10 -2 gem -2 (the line core forming region). The 
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Figure 6. Intensity formation diagram for the two Mg II lines around 279.88 nm. Cases (a) and (b) correspond to two different columns from the flaring 
simulation. The format description is the same as for Figure 1. Note the different height and temperature scales between (a) and (b). 


columns with emission lines appear clearly clustered in re¬ 
gions with a large temperature difference and high electron 
density (the distribution of n e for emission columns peaks at 
around 7 • 10 17 m~ 3 ). 

4.2. Emission as a quantitative diagnostic 

The presence of emission in the triplet lines can be an 
important indicator that the lower chromosphere has been 
heated. In addition, it can also be used to quantify the tem¬ 
perature increase that leads to emission. 

We find that the Eddington-Barbier approximation holds for 
the triplet lines. Under LTE conditions one would expect the 
intensity to follow the temperature variations. However, de¬ 
spite some coupling to the local temperature (see discussion 
above), the source function does not completely follow the 
quick temperature increases that give rise to emission. Still, 
we find that one can nevertheless derive a quantitative esti¬ 
mate of the temperature increase from the line intensity, and 
demonstrate it in Figure 5. We plot an observable from the 
two blended triplet lines, the / co re /7wing intensity ratio between 
the line core maximum (in the 279.866 — 279.893 nm inter¬ 
val) and the line wing (here taken at 279.932 nm) against 
AT, a temperature difference given by different quantities. 
The different AT are differences between the line core and 
wing forming regions. A7V = T (Smax) — T m i n gives the phys¬ 
ical temperature difference in the line forming region, while 
ATs = SVnax — T m i n is a proxy for the temperature difference as 
“measured” by the source function; the discrepancy between 
the two is a measure of how the source function departs from 
the local temperature. 

As seen in Figure 5, both AT are correlated with the ob¬ 
servable / C ore//wing- When plotted against ATs, the rela¬ 
tion with log (/ core // w ing) is nearly linear - a consequence of 
the Eddington-Barbier approximation: the source function 
and intensity are closely correlated. When plotted against 


ATs, one can still mostly recover the linear relation with 
log (/ C ore/Aving)> but there is an added offset and increased 
scatter. For clarity we show only the results for the quiet Sun 
simulation in Figure 5, but the results for the flaring simula¬ 
tion are essentially the same. Within the uncertainties, this 
enables one to use /core Aving to directly quantify the localized 
heating in the lower chromosphere when the triplet lines are 
in emission, making it a powerful diagnostic. 

In Figure 6 we show example line formation diagrams for 
the flaring simulation. Case (a) is an extreme example with 
strong velocity and temperature gradients alternating in rapid 
succession in height. The resulting triplet lines are formed 
over a very large height range and show a strong emission 
profile. There are two large peaks of the source function at 
around z = 4, 6 Mm, which roughly follow large temperature 
increases. Near the line core the z = 6 Mm peak is domi¬ 
nant, making the line profile nearly single peaked. Case (b) 
depicts a very different scenario: here the line is formed over 
a shallow height range (z < 1.3 Mm), where the velocities are 
relatively small. In this case there are two temperature bumps 
at z ~ 0.4, 0.7 Mm, where the densities are high enough that 
the source function follows. This causes a triple-peaked pro¬ 
file. There is another, much larger, temperature increase at 
z ~ 1.3 Mm and corresponding bump in the source function. 
However, at these higher layers there is very little opacity at 
these wavelengths, and correspondingly their contribution to 
the line profile is negligible. (At z = 1.33 Mm the column 
mass in this location is 6 • 10 -5 g cm -2 , just outside the range 
where the lines are sensitive.) 

As evidenced by the line profiles in Figures 1 and 6, the 
shape of the emission lines also provides important informa¬ 
tion about the underlying physical conditions. Double, triple, 
or multi-peaked profiles are mostly caused by a temperature 
profile with several rapid variations in height. Single-peaked 
lines tend to come from a dominant temperature increase, typ- 
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Figure 7. Examples of triplet emission in IRIS observations. Left panels: images from spectral rasters, taken at fixed wavelengths in the 279.88 nm lines (for 
each panel, the wavelengths were taken where the peak emission occurs at the spatial position of the red crosses). Rasters were taken in different regions and 
observing modes (see text for details), including active region (AR), near a sunspot, in the quiet Sun (QS), and in quiet Sun near the limb (QS Limb). Each 
image is oriented in solar (x, y) coordinates, and has a field-of-view of 2274 x 120"; the coordinates of the bottom left corners are shown below the images. 
Right panels: Individual spectra from selected points in the rasters on the left, labelled accordingly. Red curves correspond to the red crosses on the images, cyan 
curves to the cyan plus signs (only for AR and Sunspot). The black dashed curves depict the spatially averaged spectra over the windows displayed on the left. 


ically in higher layers. If a single temperature increase occurs 
deep in the atmosphere, at high densities, the emission man¬ 
ifests itself in the far wings of the triplet lines, with the line 
core a normal absorption profile (when this happens, the h&k 
lines are also much wider than normal). In the quiet Sun sim¬ 
ulation, there are in fact twice as many columns with far wing 
emission profiles than columns with emission in the line core. 

5. OBSERVING EMISSION WITH IRIS 

The Mg II triplet lines are routinely observed with the IRIS 
spectrograph. In Figure 7 we show some IRIS raster images 
and individual spectra of locations with emission in the triplet 
lines. The observations comprise active regions (AR, one of 
them including a sunspot) and quiet Sun (both at disk-center 
and at the limb). Details about the observations are summa¬ 
rized in Table 1 . The observations were obtained by scanning 
a region of the Sun with a moving slit, with a step size of 
0735 (“dense” rasters), except the quiet Sun dataset, which 
has a step size of 1". We made use of IRIS reduced and cali¬ 
brated level 2 data (see De Pontieu et al. 2014, for details on 
the reduction procedure). 

From the observations we find that emission in the triplet 
lines is rare. They are most easily seen in emission in flares 
and other high-energy phenomena where there is heating in 
the lower chromosphere. Vissers et al. (2015) find that the 
triplet lines are often in emission in Ellerman bombs (Eller- 
man 1917), and they are also seen in emission in some of the 
explosive events reported by Peter et al. (2014) and Schmit 
et al. (2014). 


Table 1 

Observational data sets. 


Target 

Starting time 

Center coord. 

Exp. time 

Active Region 

2014-06-15T08:29 

(419"8, 280"4) 

4s 

Sunspot 

2014-07-03T10:10 

(—302"6, -225"0) 

4s 

Quiet Sun 

2013-09-18T07:39 

(34"9, 68"0) 

4s 

Quiet Sun limb 

2014-07-10T00:10 

(538"3, -770"4) 

30 s 


The raster images in Figure 7 are taken at wavelengths close 
to the center of the lines at 279.88 nm, meaning that bright 
areas are locations of increased intensity in the lines (in the 
AR and sunspot images many of these are indeed in emission, 
but not in the quiet Sun images). In the AR images one sees a 
collection of several bright dots with a round shape where the 
lines are in emission - these generally occur in the vicinity 
of sunspots, and it is possible that they are related to Ellerman 
bombs. But as shown in the sunspot panel, there is also strong 
triplet emission (and strong h&k emission) in the light bridge. 
In some of these AR locations the h&k lines have an intensity 
very close to the average, while the triplet lines are strongly 
enhanced (see the extreme example in the sunspot spectrum 
where the lines at 279.8 nm are stronger than the k line). Such 
scenarios could indicate an abrupt temperature rise only in the 
lower chromosphere, not felt by the h&k line centers. Other 
locations show both the triplet lines and the h&k lines strongly 
enhanced, which could be caused by a temperature increase 
throughout a wider range of the chromosphere (e.g. in flares). 
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In the quiet Sun, emission in the triplet lines is seldom 
found. From several datasets investigated, it was found only 
in very few of them. And even when found, it is rarely of the 
same magnitude as seen in ARs. Most often, this quiet Sun 
emission is on the far wings of the line, with the central part 
of the line resembling its typical absorption profile. This is 
consistent with what we found in the quiet Sun simulation, 
and suggests that in these locations the heating is limited to 
a deep area near the temperature minimum. Another piece of 
corroborating evidence for this scenario is that the h&k lines 
are also noticeably wider than the average, an indication that 
the ki/hifeatures are being formed lower down in the atmo¬ 
sphere as a result of the chromospheric temperature increase 
taking place deeper than usual. This is the case with the quiet 
Sun profiles that we show - the triplet lines are in emission 
only in the wings, and the h&k are wider than normal. These 
locations also show an enhanced photospheric temperature, as 
evidenced by the high local continuum. 

6. CONCLUSIONS 

We have studied the formation of the Mg II triplet of lines 
that lie near the h&k lines. The lower levels of these subordi¬ 
nate lines are the upper levels of the h&k lines. To understand 
their formation we use a forward modeling approach, making 
use of realistic 3D radiative MHD simulations and compar¬ 
ing the predicted spectra with the physical quantities from the 
simulation. 

In the quiet Sun, we find that the lines are formed just above 
the temperature minimum, at heights around 0.6 — 1.2 Mm 
above T 500 = 1, or at column masses down to 10 -4 gem -2 . 
The lines can be used as velocity indicators for those layers, 
but the line pair at 279.88 nm is problematic because of the 
overlapping lines. The source functions of these lines typi¬ 
cally decouple from the local temperature before the chromo¬ 
spheric temperature rise, and therefore the source function de¬ 
creasing with height gives rise to absorption lines, by far the 
most common scenario. However, under particular circum¬ 
stances the lines go into emission, and this can be a powerful 
diagnostic. 

When a large temperature increase is present in the lower 
chromosphere the source functions can follow this increase, 
leading to emission lines. From the simulations we find that 
this typically happens when there is a temperature increase 
of more than 1500 K in layers with column masses from 1 
to 10 -3 gem -2 , and an electron density above 10 17 m -3 . In 
addition, one can use the ratio of the emission peak to the lo¬ 
cal continuum of the lines to derive a rough estimation of the 
temperature difference that caused that same emission. This 
holds true for both quiet Sun and more violent flaring sim¬ 
ulations, and means that the lines can be used to diagnose 
steep temperature increases in the lower chromosphere, a new 
type of diagnostic complementary to those of the h&k lines 
(Leenaarts et al. 2013a,b; Pereira et al. 2013). The shape of 
the emission line also provides information about the under¬ 
lying physical quantities: in cases where the heating occurs 
deeper down and is covered by cooler material, the emission 
in the triplet lines occurs predominantly in their far wings, 
with the central part of the line being like a typical absorption 
line (under these circumstances the h&k lines are also wider). 
When the heating occurs higher in the column range to which 
the lines are sensitive (10 -2 — 5 • 10 -4 gem -2 ), the emission 
takes place in the line core. 

The Mg II triplet lines are routinely observed by IRIS, and 


we find several example observations that confirm the scenar¬ 
ios seen in the synthetic spectra. As in the quiet Sun sim¬ 
ulation, emission is very rare in the quiet Sun, and when it 
happens it tends to be in the far wings of the lines, suggesting 
heating occurring deeper in the chromosphere. The lines are 
more easily observed in emission in flares, active regions, and 
in particular near sunspots and features like Ellerman bombs, 
as has already been reported. In such locations the lines can 
be strongly enhanced, in some extreme cases even stronger 
than the h&k lines. With the help from the simulations, one 
can now understand better some of these phenomena and use 
the diagnostics from these lines to trace instances of strong 
heating in the lower chromosphere. 
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